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§ ' Abstract 

Gravitational properties of a hedge-hog type topological defect in two extra dimensions are con- 

O ■ 

. sidered in General Relativity employing a vector as the order parameter. The developed macro- 

scopic theory of phase transitions with spontaneous symmetry breaking is applied to the analysis 
Q ' of possible "thick" brane structures. The previous considerations were done using the order pa- 

^ ■ rameter in the form of a multiplet in a target space of scalar fields. The difference of these two 

approaches is analyzed and demonstrated in detail. 

There are two different symmetries of regular solutions of Einstein equations for a hedgehog type 
vector order parameter. Both solutions are analyzed in parallel analytically and numerically. Reg- 
ular configurations in cases of vector order parameter have one more free parameter in comparison 
with the scalar multiplet solutions. 

It is shown that the existence of a negative cosmological constant is sufficient for the spontaneous 
symmetry breaking of the initially plain bulk. 

Regular configurations have a growing gravitational potential and are able to trap the matter 
to the brane. Among others there are solutions with the gravitational potential having several 
points of minimum. Identical in the uniform bulk spin-less particles, being trapped within separate 
points of minimum, acquire different masses and appear to an observer within the brane as different 
particles with integer spins. 
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I. INTRODUCTION 



The theories of brane world and multidimensional gravity are widely discussed in the 
literature. They continue numerous attempts to find the origin of enormous hierarchy of 
energy and mass scales observed in nature, to explain the dark matter and dark energy 
effects, and other long-standing problems in physics. One can find a great variety of separate 
brane-world models in the literature: thin and thick branes in five and more dimensions, 
fiat or curved branes and bulk, different mechanisms of brane formation, etc. 

From my point of view the most natural approach is to consider the appearance of a 
distinguished hyper-surface in the space-time manifold as a result of a phase transition 
with spontaneous symmetry breaking in the early phase of the Universe formation. The 
macroscopic Landau theory of phase transitions with spontaneous symmetry breaking allows 
to consider the brane world concept self-consistently, even without the knowledge of the 
nature of physical vacuum. It is reasonable to regard the brane world as a topological 
defect, which inevitably appears as a result of a phase transition with spontaneous symmetry 
breaking. 

The properties of topological defects (strings, monopoles, ...) are generally described 
in general relativity with the aid of a multiplet of scalar fields forming a hedgehog-type 
configuration in extra dimensions (see l|] and references there in). The scalar multiplet plays 
the role of the order parameter. The hedgehog-type multiplet configuration is proportional 
to a unit vector in the Euclidean target space of scalar fields. Though this model is self- 
consistent, it is not the direct way for generalization of a plane monopole to the curved 
space-time. 

In a fiat space-time there is no difference between a vector and a hedgehog-type multiplet 
of scalar fields. On the contrary, in a curved space-time scalar multiplets and real vectors are 
transformed differently. For this reason in general relativity the two approaches (a multiplet 
of scalar fields and a vector order parameter) give different results which are worth to be 
compared. It is the subject of this paper. A priory it seems more difficult to deal with a 
vector order parameter, and, probably, it is the reason why I couldn't find in the literature 
any papers considering phase transitions with a hedgehog-type vector order parameter in 
general relativity. 

The spontaneous symmetry breaking with a hedgehog-type vector order parameter in ap- 
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plication to the brane world with two extra dimensions is considered in this paper. General 
formulae (section [IT]) are applied to the analytical (section and numerical fllVp anal- 
ysis of gravitational properties of global strings in two extra dimensions. The results are 
summarized in section IVl 



II. GENERAL APPROACH 



A. Lagrangian 

The order parameter enters the Lagrangian via scalar bilinear combinations of its covari- 
ant derivatives and via a scalar potential V allowing the spontaneous symmetry breaking. If 
0/ is a vector order parameter, then V should be a function of the scalar (p^4>K = g^^4>i4>K, 
and a bilinear combination of the derivatives is a tensor 

SiKLM = 4>I;k4>L;M- 

Index is used as usual for covariant derivatives. There are three ways to simplify Sjklm 
into scalars, so the most general form of the scalar S, formed via contractions of Siklm, is 

s = A (0^)' + S0f^0f + (1) 

where A, B, and C are arbitrary constants. Different topological defects can be classified by 
these parameters. In curved space-time the scalar S depends not only on the derivatives of 
the order parameter, but also on the derivatives of the metric tensor. This is the principle 
difference between a vector and a multiplet of scalar fields. 

The general form of the Lagrangian determining gravitational properties of topological 
defects with a vector order parameter is 

, IK dgnA _ T , T 



where 



= -4 {4>!'k) ' + B't'iK't'f + C4>'k4'" - V (<!>'' 4>k) ■ (3) 

Lg is the Lagrangian of the gravitational field, R is the scalar curvature of space-time, 
is the (multidimensional) gravitational constant, and is the Lagrangian of a topological 



defect. Covariant derivation 

d(f)p I LA f 9gAM , dgAp dgMP 



'JP:M 



dx^'^ 2^ \ dxP dx^^ dx^ '^^ ^ ^ 



and razing of indexes (p^ = g^^(j)i contain g^^ and , and for this reason it is convenient 
to consider the Lagrangian as a function of 0/, g^^ , and 



B. Energy-momentum tensor 

Varying the Lagrangian dS} with respect to 5g^^ and having in mind that 

^QiK = -QKuQiN^g^''^ (5) 
we get the following expression for the energy- momentum tensor ^: 



rp _ 2 

J-IK — -JT 



dg"^ + 9QK9PIg^ I I 



(6) 



In the case of the vector order parameter the symmetry breaking potential V [(^^(Pk) = 
V (^g^^cpjcpx) also takes part in the variation with respect to Sg^^^ . 

In what follows we apply this approach to the particular topological defect, namely the 
global string in two extra dimensions. 

III. GLOBAL STRING IN EXTRA DIMENSIONS 

In my previous papers with Bronnikov (see [l| and references there in) we considered 
global monopoles and strings as topological defects with the order parameter in the form of 
a hedge-hock type multiplet of scalar fields in some fiat target space. The aim of this paper 
is to describe these defects using vector order parameter and compare the results. 

A. Metric 

The direction of the vector specifies one coordinate, and in the most simple case the 
system is uniform and isotropic with respect to all other coordinates. In our recent paper 
[ij we presented the detailed properties of global strings in two extra dimensions. For this 
reason I consider below a topological defect in the space-time with two extra dimensions. 
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The order parameter is a space-like vector {g^^(l)i(f)K < 0) directed normally from the brane 
hypersurface and depending on the only one specific coordinate, namely - the distance from 
the brane. The whole {D — dQ-\- 2) -dimensional space-time has the structure M'^^x R^x 
and the metric 



ds" = e^^Wry^.dx'^fix'^ - {df + e^^W V) , (7) 

where 77^1, = diag (1,— 1,...,— 1) is the do-dimensional Minkovsky brane metric (do > 1, we 
set do = 4 in computations) , and ip is the angular cylindrical coordinate in extra dimensions. 
7 and j3 are functions of the distinguished extradimensional coordinate / - the distance from 
the center, i.e. from the brane. e^^'^ = r [1) is the circular radius. Greek indices /x, .. 
correspond to do-dimensional space-time on the brane, and /, i^, ... - to all D = do -|- 2 
coordinates. The metric tensor qik is diagonal, and its nonzero components are denoted as 
follows: 



giK 



=27 



/ = X = 0, 



=27 



-e"', Q < I = K < do, 
-1, I = K = do, 
I ^K^^p. 

The curvature of the metric on brane due to the matter is supposed to be much smaller 
than the curvature of the bulk caused by the brane formation. 



B. Regularity conditions 

If the influence of matter on brane is neglected, then there is no physical reason for 
singularities, and the self consistent structure of a topological defect should be regular. A 
necessary condition of regularity is finiteness of all invariants of the Riemann tensor of 
curvature. The nonzero components of the Riemann tensor are 



T3AB 
^ CD 



= < 



-i^ {^M - ^D^c) . A,B,C,D<do, 
-/?Y, A^C^if, B,D<do, 

+ A = C = do, B,D<do 

- iP" + /3'2) , A = C = do, B = D = ip. 



(8) 
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Here prime denotes d/dl. One of the invariants of the Riemann tensor is the Kretchmann 
scalar K = R^^^R'^^b ; which is the sum of all nonzero R^^qd squared. Hence all the 
nonzero components of the Riemann tensor, and namely 

y, + 0i\ f3" + f3'' (9) 

must be finite. The center r = is a singular point of the cylindrical coordinate system. 
The absence of curvature singularity in the center follows from the finiteness of the last term 
in ^ . Let 

(3" + (3''^ = c < oo at / = 0. (10) 
Integrating (|TOl) in the vicinity of the center we have 

(5' = \ + \cl + [l') . (11) 

Relation (ITTll ensures the correct (= 27r) circumference-to-radius ratio, or, equivalently, 
dr"^ = dP at I 0. Finiteness of f3''y' at / = is fulfilled if 

y = (/) or smaller at / 0. (12) 



C. Vector order parameter 

Our aim is to consider the order parameter as a vector in extra dimensions directed 
normally from the Minkovsky hypersurface. In the cylindrical coordinate system of extra 
dimensions the only nonzero component of the vector order parameter is 

(f)do = 0. (13) 

In the space-time with the metric ([7]) the covariant derivative 

<Pi;K = - lSiK9''9'n<l^ (14) 

is a symmetric tensor: 0/ -k = 4>k ;/• For this reason (t)^.K4>i^ = 4'^-k4>^i 5 without loss 
of generality we set C = 0. The Lagrangian ^ takes the form 

La = A (0' + \<PY.K9''''g'KKY + B + 102 [g^^g'^^f^j - V . (15) 

It contains only two arbitrary constants A and B. In (!T5!) we set g'^'^'^^ = — 1 in accordance 
with ([7]) . However one should keep in mind that (fTSjl cannot be used in ([6]) . To derive 
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the energy- momentum tensor ([6]) one should use the Lagrangian ([3]) , and set g'^°'^° = —1, 
^gdodoy _ g g^f^gj^ differentiation. Nevertheless, the field equation can be derived using f[T5|) 
in the general formula 

' ^ y \ _ ^ ^ 0. (16) 



In the space-time with metric ([7]) the sums in f[T^ are 

^n = ^Ei.(^?^^^7W)" = rfo7'" + r, n=l,2,... . (17) 

and the determinant of the metric tensor is 

^ = (-1)^-1 e2('^°^+^). (18) 

We consider below the two CclSGS, CcLSG A (A 7^ 0, -B = 0) and case B (S 7^ 0, A = 0), in 
comparison, but separately. 

D. Field equations 



Substituting ( 1151) into (1161) we get the field equations in the case of vector order parameter. 
Case A 



In case A = i, i? = we get 



h' + ^10)' + = 0- (19) 



Case B 

In the case I? = ^ yl = the substitution results in 



dV 

f + Si(t>' -S2(t>+^ = Q. (20) 

S*! and 5*2 are defined in (ITTj) . 
Case of scalar multiplet 

n 

In the case of the multiplet of scalar fields we had jll|: 

0" + S,<j>' - 0e-2^ + ^ = 0. (21) 

0(p 

Unlike (I2T!) . the field equations (fT9l) and (!20|) do not depend directly on /? (and thus on 
the circular radius r = ln/3). The case A field equation (|T9l) includes second derivatives of 
the metric tensor. In the flat space-time 7' = 0, P' = j, P" = —j2, e'"^^ = and all three 
field equations ([19]), (ED]), and (EI]) reduce to 



ill 



+ y0'- ^0+1^ = 0. (22) 
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E. Energy-momentum tensor 



The energy-momentum tensor ([6]) inevitably contains second derivatives. However, in 
case A the second derivatives can be excluded with the aid of the field equation f[T^ . In 
case B the second derivatives cannot be excluded from the energy-momentum tensor. The 
final result of rather wearing derivations is: 

Case A 

Tf = W + S,<pf + 5fV + {5pf^ - 5f) f (23) 



Case B 



rpK 

J^do _ 
do 

rpK 
^I>do 



5f 



(24) 



K 



:^/?')' 0' + /5' (0Y - + 1^20^) + ^' 



Unlike the scalar multiplet case, the energy-momentum tensor (l23l) contains not only the 
potential but also its derivative 

Correctness of (j23|) and (!24l) is checked by the derivation of the covariant divergence Tf.j^ 
(actually .^). Again, with the aid of field equations ([T9l) and fl20l) we confirm that 



■do 



0. 



F. Einstein equations 



The same way as in [i] we use the Einstein equations in the form 



where Rf is the Ricci tensor. 



K 



5f (7" + y^i), Kdo 
51{S[ + S,), I = do 



and 



= - i^f T. 
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In both cases A and B the sets of Einstein equations consist of three first order equations 
with respect to 7', and [5^]: 
Case A 



7" + ^17' 
S[ + S2- 



2 2]/ _|_ 1 dV ^ 

do do d(f> 



K 



do 



f3" + Si(3' = 



do 



(25) 
(26) 
(27) 



Case B 



S'l + 5*2 



do 



do 



do 
do 



(28) 
(29) 
(30) 



Metric functions 7 and /3 do not enter the equations (j25l — I 30p directly, only via the 
derivatives. In the case of a scalar multiplet order parameter, see eq.(14 — 16) in /3 
enters the Einstein equations directly, and the system of equations is of the fourth order. 
It is also worth to mention that the right hand sides of first and third equations coincide 
in case A, and have a similar structure in case B. This symmetry reflects the fact that the 
coordinates with I < and I = ip are cyclic variables. 

The field equations ( |T9l) and ( l20l) are not independent. They are consequences of the 
Einstein equations due to the Bianchi identities. 

G. First integrals 



Excluding the second derivatives 7" and [3" in the sets (1251- 1 27p and (1281- 1 30p , we get 
the relations 
Case A 

51 -82 = -n" [(0' + Si0)' + 2V\ , (31) 

Case B 

52 (1 - /t^') = + 2«:2y + Si, (32) 
which can be considered as first integrals of the systems (1251 — 1 271 ) and (1281— 1 301) . 
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H. Further simplification 

The equations (125!) and (127|) , as well as the equations (l28i) and (130|) , have similar struc- 
tures. Extracting in both cases one from the other, we get the equations 
Case A 

(7 - /?')' + (7 - ^1 = 0, (33) 

Case B 

[(y - /?') (1 - + (V - (1 - ^^') S, = 0. (34) 

fl33|) can be used instead of one of the equations fl25l) or fl27j) and, respectively, fl34l) - instead 
of (EHD or (EDI) . 

With the aid of the relations (13T]) . ( l33il and ( 132I) .( 1M|1 the complete sets of equations can 
be reduced to more simple forms. Introducing new functions 

U = -f'-p', Z = 4>' + Si^, = (/.', (35) 

and having in mind that functions jS', 7', and their combination S2 = ^07'^ + P''^ (flTj) are 
expressed via U and as follows: 

, U + Si ^, Si-d^U ^ doU^ + Sf 
rfo + 1 rfo + 1 rfo + 1 



we get the sets of four first order equations, solved against the derivatives: 
Case A 

U' = -USi 
(j)' = Z- Si(f) 



Case B 



[U (1 - /€202)]' + U{1- fi:^') Si=0 



2(l+rfo) ,^2 



do 



= 



(p' = 1p 



+ =0 



(36) 



(37) 



The sets (136|1 and (l37|l are most convenient for both analytical and numerical analysis. 
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I. General analysis of regular solutions 



The sets ( 136|) and (1371) are invariant against adding arbitrary constants to 7 and /3. 
Without loss of generality we can set 

7(0) = 0. (38) 

Requirement of regularity in the center dictates the condition fllip , and, if we do not consider 
configurations with angle deficit (or surplus), we have 

r = = / at / ^ 0. (39) 

Integrating first equations in the sets (15B11 and (P7|l with boundary conditions (13811 39p we 
get: 

Case A 

7'-/5' = -e-(''«''+^), (40) 

Case B 



V-/'' = -^-i-i. (41) 



-(d07+/3) 

From (HOj) and (14T|) it follows that in both cases (3' > 7' everywhere. Within the regions of 
regularity e"^'^"'''"'"^^ at / ^ 00, and 7^ coincides with 

The denominator in (|4T|) is the origin of an important difference between the cases A and 
B. In the symmetric (unbroken) state the order parameter = 0, and 0' = everywhere. 
In the state with broken symmetry 0' 7^ 0, but still = at a point of maximum of 
the symmetry breaking potential V(0), which is supposed to coincide with the center / = 
0. Requirement of regularity is fulfilled only if the denominator in fHTj) remains positive 
everywhere. It means that in case B the order parameter 0(/) can never exceed : 

< 1, case B. (42) 

Case A is free from such a restriction. 

n 

Recall that topological defects, formed as multiplets of scalar fields pj] , are of three types. 
Integral curves can terminate with: 

A) infinite circular radius r (/) at / —>■ 00; 

B) finite circular radius r^o = r (00) = const < 00; 

C) second center r = at some finite I = Ic- 
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In the vector order parameter case the situation is different. Equations (j^Oj) and (HTj) 
allow to prove that regular configurations in both cases A and B can not terminate neither 
with a finite value of circular radius Too at / — oo, nor in the second center. 

Suppose, for instance in case A, that r^o = const < oo. Then /3'(oo) = 0, and fHUl) reduces 
to 7' = — Lg-rfo7 at / — s> 00. After integration we get 

r 

Iq is a constant of integration. The l.h.s. is obviously positive, while the r.h.s. becomes 
negative and infinitely large at / 00. Thus Too = const < 00 is impossible. 

The second center is also impossible. In the vicinity of the second center the l.h.s. of 
f l40|) becomes large positive due to — and the r.h.s. remains negative. 

We come to the conclusion that regular configurations of topological defects with the 
vector order parameter start at the center / = and always terminate at / ^ 00 with 
infinitely growing circular radius r (/) 00. 

It follows from the requirement of regularity ( |T2|) that 7' = 70'/ at / — 0. From the first 
integrals (1311) and (l32l) we find the relations between 7q , 0q, and Vq : 

Case A 

7;' = -^(20;^ + ^o), (43) 
"0 

Case B 

where Vq is the value of the symmetry breaking potential at the center / = 0. In both cases A 
and B, as well as in the case of scalar multiplet, the value 0g = 0' (0) is not restricted by the 
equations. The difference is that in the scalar multiplet case 0o becomes fixed unanimously 
by the requirement of regularity, while in both cases of vector order parameter 0g remains 

a free parameter. 



J. Asymptotic behavior 

Condition of regularity requires that 7' is finite everywhere. Within the area of regularity 
it tends to a fixed finite value 7^ at / 00. As soon as r (Z) ^ 00 at / ^ 00, we see from 
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and gl]) that 7' - ^ in both cases A and B. Thus /3' (00) = 7^^. The field (/) 
also tends to its finite value (poo = (p (00) . 
Case A 

From the field equation f[T^ it follows that ^ ^ at / ^ 00, i.e. in the case A reg- 
ular configurations terminate at an extremum of the potential V{(f)). Let Voo = V {(poo) , 
V {(poo) = 0. From the first integral (13T!) we find the limiting value 7^ : 



^1 {do + l)[do + {do + l)f^^^'^' ^^^^ 



ooj 



A necessary condition of existence of regular configurations of topological defects with 
the vector order parameter is Voo < 0. 
Case B 

At / — > cx) the second equation in (1371) reduces to 

^^-^)-^^'^''^-- 

We see that in this case also 

V;o = ^(0oo)<O. (46) 
The fourth equation in fl371) reduces to 

_ C.2 0OO , dVoo _ 
'^do + l^ d(P 

and so 



Here Voo and are the values of the symmetry breaking potential and its derivative at 
/ — s> 00. Unlike the case A, in case B regular configurations terminate at a slope of the 
potential, and not at a point of extremum. Excluding Sf we get the equation determining 

(poo = (p{oo) : 

n'Vo. = -{do + n'(Pl)^. (48) 
In view of 5*100 = {do + 1) 7^, we find 

y=j ' '-Y^ (49) 

V K + l)0oo d(p- ^ ^ 
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To find the behavior of (/) and 5*1 (/) at large distances from the center we have to 
linearize the sets of equations fl36|) and fl37|) at / ^ oo : 



Case A: 

6<f)' = 6Z- {do + 1) i^6<j) - 0oo55i (50) 
6Z' = 

Here primes denote derivatives d/dl, {5S[ = d{5Si) / dl, ...) , except = |^ \<j>=(i>oo ■ Ex- 
cluding 5Z and 55*1, we get the second order linear homogeneous equation for 50 : 

2^2 IT/ \ V" 

5f + (rfo + 1) i^H' + , u T\^^>2 ^^ = 0- 

do [do + 1) l'^ 

In case A the solution terminates at an extremum of V{(j)). If the extremum of the 
potential is minimum, (V"^ > 0) , the nontrivial solution vanishes at / ^ oo : 

6(j) = Ae^+^ + Be^-\ (51) 

where A and B are constants of integration, and both eigenvalues 



/OO 



are either negative, or have negative real parts. Absence of growing solutions is the reason 
why 0Q remains a free parameter in the vector order parameter case. 

In case A the asymptotic behavior of the field (/) far from the center is determined by 
two constant parameters of the symmetry breaking potential near its extremum, namely V^o 
and V^. If the extremum is minimum, > 0, then the expression under the root can be 
both positive and negative. So (Z) can tend to 0oo either smoothly, or with oscillations. In 
the space of physical parameters the boundary between smooth and oscillating solutions is 
determined by the relation 



/oo 



do {do + 

Oscillating behavior of the field (/) induces oscillations of (3' and 7'. If 7' changes sign, 
then 7 (Z) can have minimums. Remind, that in the weak gravitation limit 7 acts as a 
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gravitational potential, so the matter can be trapped near the minimums of 7 (/). Trapping 
of matter to the brane is considered below in sections IIVBI and IIV CI in more detail. 

Usually = is a maximum of the potential V (0) . It is also an extremum, dV/ 90 = at 
= 0. Regular configurations, starting from the center I = with (0) = 0, can terminate 
at / — > 00 with 0OO = as well. In this case = V" (0) < 0, and the hnear set (|5UI) reduces 
to the following asymptotic equation for (/) : 

cp" + (do + i)y^<p'-\v:i\<p = o. 

Its general solution is a linear combination of vanishing and growing functions: 



Requirement of regularity demands to exclude the growing solutions from the consideration. 
It can be done at the expense of 0o. In case A regular solutions terminating at a maximum 
of the potential can exist only at a discrete sequence of values of 0q. 

Case B: 

Linearizing the set (1371) we get the following equation for 50 : 

V <90^ (do + 1) (do + x^0So) / 
Again its nontrivial solution has the form 

50 = Ae^+^ + Be^-\ (54) 

where A and B are constants of integration, and the eigenvalues are 



, _ "-u I - ' 00 I / -"U I - --u Y'oo ^ ' 00 ' 00 /ppx 

do + >c-<Pl^)^~^- ^^^^ 

Eigenvalues (155|) do not look as transparent as ( l52i) . However, in practice both eigenvalues 
f l55p are complex with negative real parts. Actually, in view of 



4 do + J 0oo90 902 V 4 J (t)^d(t) 902 ■ 

For the "Mexican hat" potential (|60l) . that we use in computations (see section HVl below) . 
the r.h.s. is negative at least for d^ < 7. In case B the order parameter 0(/) approaches its 
limiting value 0oo with oscillations. 
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K. Boundary conditions 

The complete sets of equations determining the structure of topological defects in the 
cases A f l2^ 27| 311) and B fl28| 30| 321) of vector order parameter are of the third order with 
respect to three unknowns 7',/^', and 0. The simple solutions are determined unanimously 
by fixing the values of these three functions in any regular point. The center / = is a 
singular point of the cylindrical coordinate system. The condition (0) = fulfills for both 
symmetries (high and broken). /3' is infinite at / = 0. For this reason we set the boundary 
conditions very close to the center, but not exactly at Z = 0. 

For numerical analysis it is convenient to deal with a system of four first order equations 
solved against the derivatives: (I3B|) in case A, and (I57|) in case B. 

Case A 

The symmetry breaking potential V (0) enters the equations (l36l) only via its derivative 
If we leave only the main terms in the boundary conditions: t/ = — y, 5*1 = y at 
/ 0, then we loose the information about the absolute value of the potential. The value 

= y (0) appears in the next approximation. Using the expansion (1111) of (3' in the vicinity 
of the center and the equation (HUj) , we express c via 7q : 

c=-K-2)7o. 

To preserve the complete information about the symmetry breaking potential one has to 
write the boundary conditions at / ^ as follows 

f/ = |(rfo + l)7o^-T, ^i = |K + l)7o/ + T, = 00^, Z = 2cl>',. (56) 

The values 79 , 0o, and Vq are not independent. They are connected with each other by (l43l) . 
Case B 

In case B the potential V (0) enters the equations ( 1371) directly. The information about 
its absolute value is not lost even if we use only the main approximation in the boundary 
conditions at / ^ : 

t/(l-xV) = -y, ^1(1 + ^) =7' ^ = ^'0^^ ^ = '^o- (57) 
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L. Case A analytical solution with 



dV 



If the potential V = Vq does not depend on 0, then it actually plays the role of the 
cosmological constant A = k^Vq. The peculiarity of the vector order parameter in case A is 
that the equations (1361) loose the information about the potential completely if ^ = 0. The 
absolute value Vq is present only in the boundary conditions (!56|) . The equations (136|) with 
^ = and boundary conditions fl56|) have the following analytic solution 



U = ^P-^, S, = VCcoth(VCl], 0(/) = Mtanh^, 

sinh(v^/j ^ ^ y/C 2 

where 

C = 2 (rfo + 1) 7^' = -^^^^ (2«:V'o' + A) . (58) 

"0 

The solution is regular if C > 0, i.e. A < — 2k^0q^. For ^foo = e^"^ and r = we find 
fl-oo (0 = e^^ = I cosh ^^^^ , r{l) = ^ f cosh 



The inclination (h'r, remains arbitrary. If S'r. = this solution reduces to the one found earlier 

n n 

(see [1| and [3[) for the special case = 0. The point is that the Einstein equations with a 
negative cosmological constant have a nontrivial solution (with a nonzero order parameter) 
even without a symmetry breaking potential. 

In case A a necessary condition of regular solutions with broken symmetry is the existence 
of extremum points of V (0) , where ^ = 0. In case V = const the condition ^ = is 
fulfilled identically, and formally the order parameter can tend to any (p^o as I oo. 
The displayed above analytical solution shows that the existence of a negative cosmological 
constant is sufficient for the symmetry breaking of a uniform plain bulk. 

The special case C = +0 in flSHl) when 7q = and 



Ml— (59) 



corresponds to the plain bulk goo (/) = 1, and r (/) = /. 
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IV. NUMERICAL ANALYSIS 



A. Regular solutions in the space of parameters 



The numerical integration of the sets of equations ([ 
'Mexican hat" potential taken in the same form as in [l 




and (1571) is performed for the 




2 



V 



4 



(60) 



The dimensionless parameter e moves the "Mexican hat" up and down. It is equiva- 
lent to adding a cosmological constant. The energy of spontaneous symmetry breaking is 
characterized by rf/'y^~'^\ and 



determines, as usual, the length scale. In most cases a is associated with the core radius of a 
topological defect. The strength of gravitational field is characterized by the dimensionless 
parameter 



Without loss of generality we set in computations a = 1 and rj = 1 (i.e., measure I in units 
a, and - in units rj.) 

In both cases A and B of vector order parameter the state of broken symmetry is controlled 
by four dimensionless parameters dQ,e, F, and 0q. It is the main difference from the scalar 
multiplet case where the regular configurations with given (io, and F existed only for some 
fixed values of 0g. Now the regular configurations with given do,e, and F exist within the 
whole interval < 0o < 0g The upper boundary of the interval 0q ^^.^ is a function of 
do,e, and F. This additional parametric freedom allows to forget about the so called "fine 
tuning" of the physical parameters. 

For visual demonstration it makes sense to fix c/q = 4 and one of the three other parame- 
ters. Then the area of existence of regular solutions can be presented as a map in the plane 
of two remaining parameters. 



In case A regular solutions terminate at the points of extremum of V(0). In dimensionless 
units the potential (l60l) has three extremum points - a maximum at = 0, and two minima 



F 



2 2 
K 7] . 



(62) 



Case A 
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It 








FIG. 1: Mexican hat potential (dot line, e = —5). In case A regular solutions terminate at the 
points of extremum (black points). In case B the solutions terminate on the slopes of the potential 
(black solid lines). 
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FIG. 2: The area of regular configurations in the plane (e, </>o) for do = 4 and F = 1. The upper 
curve is the boundary of existence of regular solutions. Other curves separate the regions with 
different signs of (poo- They quickly condense to the upper curve. Below the lowest curve does 
not change sign. 

at = ±1, (black points in Fig[TJ) At the limiting values of the order parameter 

VL = 0, \C = 2, = ±1 

vL = 0, = -1, = 0. 

In case A the value of F is not restricted from above. Fig [2l shows the area of regular 
configurations in the plane {e, 0o) for do = 4 and F = 1. Depending on the values of e and 0q 

19 



1 14 




FIG. 3: Case A. Map of regular solutions in the plane (F, — e) for (p^ = zty— do = 4. The 
solid curve (j64p separates the regions of smooth (above) and oscillating (below) behavior of the 
order parameter at / — > oo. To the left of the dash curve the order parameter changes sign. 



the order parameter (/) tends to +1, 0, or — 1 as / — » oo. The sequence of curves fn {£) iu 
Fig. [2] are those where (Z) as / oo. They separate the areas with different signs of 
0OO- Below the first curve fi (e) from the bottom, where < 0q < fi (e) , the order parameter 
(j) (l) doesn't change sign. Between fi (e) < 0q < /2 (e) it changes the sign once. In the area 
/2 (e) < 00 < /s (e) it changes the sign twice, and so on. The curves /„ (e) quickly condense 
to the upper curve /oo (^) — ^ oo. (5) is the upper boundary of existence of regular 
solutions (for the particular values do = 4 and F = 1). 
The curves in Fig. [2] are those where 



t£,cio = 4,F = l)=0. 



(63) 



Similar curves can be shown for fixed 0o in the plane {e, F) . For instance, the dash line in 

—^,t, Cto 



Fig. [3] is the first one of the curves 0oo yt'o = i 
parameter tends to zero at / ^ 00. The value (p^ 



4,F 



0, where the order 



= iy— ^ fl59l) corresponds to 7q = in 
( l43ll . It is the case C = in ( |58l) . so that the symmetry breaking of the plain bulk is caused 
completely by the potential V (0) , and not by the cosmo logical constant. To the right of 
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FIG. 4: Oscillating solutions in the vicinity of the lower black point on the dash curve in Fig. [3l 
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FIG. 5: Smooth solutions in the vicinity of the upper black point on the dash curve in Fig. [3l 

the dash hne (/) does not change the sign. 

For the potential fl60|l the boundary hne fl53l) between osciUating and smooth (p (Z) is 



£b = Id — 1 



(64) 



G ' rfor 

It is presented in Fig. [3] (solid line). Below the solid line the order parameter (/) tends to 
its limiting value (p^o with damping oscillations (see Fig. H]), and above this curve - without 
oscillations, see Fig. [5l 

The curves in Fig. H] correspond to the close vicinity of the lower black point on the dash 
curve in Fig. [3l the curves in Fig. [5]- to the vicinity of the upper black point. 
Case B 

There are two major qualitative differences between cases A and B. In case B: 
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FIG. 6: Case B. The upper boundary of regular solutions found numerically. 

1. Parameter F ( |62l) . characterizing the strength of gravitational field, is restricted from 
above, as it follows from ( H2l) . In the plane of parameters (F,e) the upper boundary of 
existence of regular solutions, found numerically, is presented in Figure El 

2. Regular solutions terminate with 0oo on a slope of the potential V{(j)), and not at the 
points of extremum. See black solid parts on the dashed curve in Fig{Tl where the conditions 
(gnD and (HZD are fulfilled. 

The final value 0oo obeys the equation ( l48l) . In case of Mexican hat potential ( l60l) we 
find 

If F is small, fl^ reduces to 

0L = i + ^, r«i, (66) 

and 0OO goes to a minimum of V{(f)) in the limit F ^ 0. 

Typical map of regular solutions in the plane of parameters (0q, F) is presented in Figure 
[71 Here e = — 10, and = 4. The solid line in Figure His the upper boundary of existence of 
regular solutions. The dotted hne separates the solutions with the same sign of (below) 
from the solutions with </>(/) changing sign ones (between dotted and dashed lines). Between 
dashed and dadot lines 0(/) changes its sign twice, and so on. There also are solutions with 
changing sign more than three times, but other separating lines are not shown in the 
Figure [71 
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FIG. 7: Map of regular solutions in the plane of parameters (0g,r). Case B, e = —10, do = 4. The 
solid line is the upper boundary of existence of regular solutions. Dotted, dashed and dadot lines 
separate the regions where 4>{l) changes its sign zero, once, twice, and three times. Separating lines 
between the regions where changes sign more than three times are not shown. 




FIG. 8: Case B. F = 0.05 and e = —5. Order parameter of the same sign (solid curve), 
changing sign once (dotted curve), and twice (dashed curve) Corresponding values of 0'(O) are 1, 
2, and 3. In accordance with ([65]) (j)oc = ±1.01530. 



The order parameter 0(/) is presented in Figure [8] for F = 0.05, e = —5, and 0'(O) taken 
from the different regions in Figure [7l Solid, dotted, and dashed curves correspond to the 
order parameter of the same sign, changing sign once, and twice. The values of (f)'{0) are 1, 
2, and 3, respectively. 
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B. Neutral quantum particle in the space-time with metric ([7]) 



A neutral spinless quantum particle is described by a scalar wave function x with the 
Lagrangian 

Lx = \9^\:bX,a - \nlx'X- (67) 
In the uniform bulk (while the symmetry is not broken) it is a free particle in the D- 
dimensional space-time with mass mo and spin zero. In the broken symmetry space-time 
with metric ([7j) it satisfies the Klein-Gordon equation 

^ (v^/'^X.a) ^ + mgx = 0. (68) 



The metric ([7]) depends on the only one coordinate x'^" = /. So the momenta, conjugate 
to all other coordinates, are quantum numbers. The wave function in a quantum state is 

X (a;^) = ^ (0 exp {-ip^x^ + in^p) , (69) 

where = (-E, p) is the cio-momentum within the brane, and n is the integer angular mo- 
mentum conjugate to the circular extradimensional coordinate y?. X (/) satisfies the equation 

111 

X" + SiX' + (p2e-27 _ n\-^^ - m^) X = 0. (70) 

The eigenvalues of = — compose the spectrum of squared masses, as observed 
within the brane. Quantum number n is the integer proper angular momentum of the 
particle. From the point of view of the observer in the brane it is the internal momentum, 
identical to the spin of the particle. 

The equation (17U]1 takes the form of the Schrodinger equation 

y..+ [p'-yg{x)]y = (71) 

after the substitution 



dl = e^dx, X(l) = y (x) / VTM, / (^) = exp {-\ [{do - 1) 7 + • 

The gravitational potential 

V,W^e-(e-V + ,„S)4^i^(jij|) (72) 
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FIG. 9: Case A. A solution with oscillating order parameter (/>(/). Here do = 4, e = —2, T = 10, 

0'(o) = 

determines the trapping properties of particles to the brane. The trapping is insured by the 
exponentially growing warp factor e^^: in view of (H5l) and (H9l) in both cases A and B 7^ 
are positive constants. In terms of U, Si, and (135!) the dependence of the gravitational 
potential (I72|) on the distance / is 

^ °^ 4 (4 + 1) 2 V -^^"^ + 1 ; 

(73) 



C. Oscillations 



The oscillations of the order parameter 0(/) give rise to the oscillations of the gravitational 
potential ( 172|) . see Figures iQllTOl 

In terms of fl64|) the eigenvalues ([52]) are 



A 



± 



8(G+ 1) 



1±a/1 + ^(G + 1)' 



(74) 



The less is |£:| the more oscillations display themselves. In the limiting cases of small and 
large F the frequencies of oscillations 
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FIG. 10: Case A. Gravitational potential (/) (I73p for the same set of the parameters as in FiglHl 
do = 4, e = —2, r = 10, = \J The initial mass of a test particle is set mo = 0. The solid 
curve corresponds to the angular momentum n = 0, and the dashed one - to n = 1. 



J v^, r ^ 0, 
(i + i)r, r^oo 

do not depend on e as / ^ oo. At |£| ~ 1 and F ^ 1 the gravitational potential has many 
points of minimum, see Fig. dUl 

The length scale a ( 1611) remains an arbitrary parameter of the theory. The physical 
interpretation is different in the limiting cases of large and small a. If a is extremely large, 
each minimum of the potential V{1) (!73l) forms its own brane. If the potential barrier is 
high, the branes are separated from one another. 

In the opposite limit, when the scale length a is extremely small, all points of minimum 
are located within one common brane, and in the spirit of Kalutza-Kline the points of 
minimum are beyond the resolution of modern devices. 

Low energy particles can be trapped by the points of minimum of the potential (!72l) . 
Identical in the bulk neutral spin-less particles, being trapped in the different minimum 
points, acquire different masses and angular momenta. If the scale length a is extremely 
small, then they appear to the observer in brane as different particles with integer spins. 



26 




FIG. 11: Case B. Gravitational potential (l) (I73p for the same sets of the parameters as for 
the curves of order parameter in FiglHJ do = 4, e = —5, F = 0.05, (j)'{0) = 1 (solid curve), 2 
(dotted curve), and 3 (dashed curve). The initial mass of a test particle mo = and the angular 
momentum n = 0. 




FIG. 12: Case B. Gravitational potential Vg (l) (j73p for a particle with mass mo = and angular 
momentum n = 1. The sets of the parameters are the same as in Figllll do = 4, e = —5, F = 0.05, 
0'(O) = 1 (solid curve), 2 (dotted curve), and 3 (dashed curve). The difference between angular 
momenta n = and n = 1 is essential only near the center. 

V. CONCLUDING REMARKS 

As a result of rather wearing derivations of the energy-momentum tensors we get more 
simple equations than in the case of scalar multiplet models. 
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TABLE I: Mutual comparison of main features of regular solutions employing the vector order 
parameter (cases A and B) with the scalar multiplet model, ny is the number of free parameters 



of the symmetry breaking potential V{(p). 








T" o n T" "i" v 

X J. L V 




Case B 


Scalar multiplet 






(J 


3 


4 






TJ T ^ -1- 1 


ny + 1 


nv 




L 111c L Lllllllg, 


1 iU 1 I/OOUj 


no need 


in some cases 




r(oc) 


oo 


oo 


oo,rm,0 




Matter trapping 


yes 


yes 


yes 




Presence of dV/d(p in equations 


yes 


no 


no 




Presence of V in equations 


no 


yes 


yes 




Derivation of Tjx 


wearing 


wearing 


easy 




Equations are 


most simple 


more simple 


simple 




Strength of gravitational field T 


arbitrary 


restricted from a 


hove arbitrary 



The main features of spontaneous symmetry braking with a hedge-hog type vector order 
parameter in comparison with the widely used previously scalar multiplet model are pre- 
sented in table [H The solutions have additional parametric freedom: in both cases A and B 
the inclination 0g is arbitrary within a whole interval of values (see Figures [2] and [6]) , and 
not only at some fixed values as in the scalar multiplet case. It means that the possibility 
of existence of the brane world is not connected with any restrictions of fine-tuning type. 
The origin of the additional parametric freedom is the order of equations, which in case of 
vector order parameter is less than in scalar multiplet models. 

All regular configurations display trapping properties. Oscillating behavior of the order 
parameter, especially in case A, gives rise to existence of several points of minimum of the 
attractive potential (1731) at different energy levels. Particles, trapped at different points of 
minimum, acquire different masses. If we assume that the length scale a (16T1) is extremely 
small, it can be a reason of the observed hierarchy of masses. Angular momentum of extra- 
dimensional motion n is also a quantum number. It appears to the observer on brane as an 
internal momentum of the particle which can be scarcely separated from its spin. 

Most elementary particles have half-integer spins. The simple case of spontaneous sym- 
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metry breaking, considered above, can not connect the origin of half-integer spins with 
extra-dimensional angular momenta. Half-integer spins in General Relativity is still a prob- 
lem. 
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